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ABSTRACT 

Chandra observations of elliptical galaxies have revealed large numbers of Low Mass X-ray Binaries (LMXBs) 
accreting at rates > 10~^Mq yr~'. One scenario which generates this M from an old stellar population is nuclear 
driven mass transfer onto a neutron star or black hole from a Roche lobe filling red giant. However, in our Galaxy, 
most of these systems accrete sporadically as transients due to a thermal instability in the accretion disk. Using 
the common criterion for disk instability (including irradiation), we find that this mode of mass transfer leads 
to transient accretion for at least 75% of the binary's life. Repeated Chandra observations of elliptical galaxies 
should reveal this population. The recuiTence times might be very long (~ 100-10,000 years depending on the 
orbital period at the onset of mass transfer), and outbursts might last for 1-100 years. Mass transfering binaries 
can also be formed in old populations via interactions in dense stellar environments, such as globular clusters 
(GCs). These tend to have shorter orbital periods and are more likely stable accretors, making them apparently a 
large fraction of the elliptical's LMXB population. 

Subject headings: accretion, accretion disks: binaries: close — galaxies: elliptical and lenticular — pulsars: 
individual: GRO J1744-28 — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Chandra observations of nearby elliptical and SO galaxies 
(e.g. Sarazin, Irwin & Bregman 2000; Blanton, Sarazin & Irwin 
2000; Sarazin, Irwin, & Bregman 2001; Finoguenov & Jones 
2001; Angelini et al. 2001; Kraft et al. 2001) often find « 20 
X-ray sources with luminosities greater than 10^^ erg sec"' and 
even more with luminosities greater than lO-'^ erg sec"'. For a 
neutron star (NS) accretor, this corresponds to accretion rates 
as large as 10"'^Mo yr"', implying a short lifetime compared 
to the age of the stellar population. One way to generate such 
rapid accretion rates in an old stellar population (neglecting for- 
mation of tight binaries in dense environments, such as globular 
clusters) is to have a Roche lobe-filling binary driven by nuclear 
expansion of a low-mass star on the red giant branch (orbital pe- 
riods are days or longer). Ritter (1999) has shown an analytic 
solution for such a Low Mass X-Ray Binary (LMXB), which 
is further simplified in ellipticals where the initial companion 
masses should all be k, 0.9Mq for an old population. 

However, such wide binaries are subject to thermal insta- 
bilities in the accretion disk (King, Kolb & Buderi 1996; see 
Dubus, Hameury & Lasota 2001 for an up-to-date discussion) 
that cause transient accretion. Using the cuiTent instability cri- 
terion, we follow the evolution of such an LMXB and find that 
over 75% of its life is spent as a transient. We thus expect that a 
large fraction of the field LMXBs in elliptical galaxies are long- 
term transients. Our ability to detect them with repeated Chan- 
dra observations depends on their outburst durations, which can 
be as long as years. 

We overview the binary evolution scenario in §2 and sum- 
marize the challenge of finding the bright, blue optical coun- 
terparts to these LMXBs in elliptical galaxies. In §3, we apply 
the current understanding of accretion disk instabilities in X-ray 
sources and estimate the outburst recurrence times. In §4, we 



compare our work to galactic systems, noting that the recently 
discovered galactic transients with outburst durations in excess 
of ten years (e.g. Wijnands 2002) might well be wide binaries. 
We close in §5 with our conclusions and a discussion of alter- 
native binary scenarios that lead to more stable mass transfer, 
especially relevant for those LMXBs found in globular clus- 
ters (Sarazin, Irwin & Bregman 2001; Angelini, Loewenstein 
&Mushotzky 2001). 

2. BINARY EVOLUTION WITH RED GIANT DONORS 

We only consider binaries with initial orbital periods > 1 day. 
If the binary is tighter (Porb < 1 day), rapid accretion may be 
initiated due to orbital angular momentum loss from magnetic 
braking, resulting in mass transfer while a low-mass star is still 
on the main sequence (see Tauris & Savonije 1999). Whenever 
there are multiple parameters to choose, we always take the 
value that maximizes the stability of the accretion disk. This 
assures that our calculated fraction of the evolutionary lifetime 
spent as a transient is a firm lower limit. 

We use Ritter's (1999) analytical solution for the LMXB evo- 
lution, which assumes that the donor star is Roche lobe filling 
on the red giant branch and expanding due to an increasing lu- 
minosity. The mass transfer is slow enough that the donor star 
always remains close to thermal equilibrium; thus Ri = R2,r, 
where R2 is the radius of the donor star and R2.R is the radius of 
its Roche lobe. For simplicity, we assume that mass and angu- 
lar momentum are conserved, which implies that at the end of 
the LMXB's evolution the NS (of mass Mi) will have a mass 
of 1.7-2. 1M0. We assume that the initial donor mass is O.9M0 
with a population I metallicity, consistent with recent determi- 
nations of the age and metallicity of the predominant stellar 
populations in these elliptical galaxies (Trager et. al. 2000). 

Once the evolution is solved, all parameters of the binary 
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can be found in terms of the orbital period at the onset of mass 
transfer, f,, which is measured in days, and the donor's mass, 
M2, which is measured in solar masses. The mass transfer rate 
within the binary is 
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where 1 AMq is the initial NS mass and 23Mq is the total mass 
of the binary. Since we will show in §3 that the disk is likely un- 
stable, we will not concern ourselves here with the implied su- 
per Eddington accretion rates at large f , should this mass trans- 
fer rate proceed onto the NS. Using equation (|l]) the evolution 
of the binary system is shown for varying P, by the solid lines 
in Figure 1. The long dashed line is the boundary where the 
hydrogen envelope of the donor has vanished, halting the evo- 
lution at the ending orbital period noted. The evolving orbital 
period is 
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shown in Figure 2 by the soUd Unes. 

For such an LMXB in an elliptical galaxy, we expect: (1) the 
orbit must be wide in order to reach such high M's, and (2) the 
reprocessing of the X-rays by the large accretion disk should 
make it a bright, very blue object. At the Eddington luminos- 
ity, an X-ray binary with Porb =10 days (much like Cygnus 
X-2) has Mv = -2 (van Paradijs & McClintock 1994) and is 
blue (B-V « Q,U-B « -1; van Paradijs & McClintock 1995). 
For NGC 4697 (at a distance modulus of « 30.35; Tonry et al. 
2001), this corresponds to an apparent V magnitude of 28.4. 

Assuming that the LMXB is not associated with a globular 
cluster, the brightest star which could confuse a search for the 
optical counterpart would be a red giant branch star just about 
to reach helium core ignition (a TRGB star). For metallicities 
[Fe/H] < -1, Mv(TRGE) w -2.5, whereas for [Fe/H] « -0.7, 
MviTRGB) « -1.5 (Saviane et al. 2000). For a slightly metal- 
rich star (such as is likely the case for these ellipticals; Trager 
et al. 2000), the peak brightness is even fainter, MviTRGB) « 
-1.0 (Garnavich et al. 1994). Hence, the optical counterpart to 
a wide LMXB will be brighter at V (and especially so in B or 
U) than any red giant branch star in an elliptical galaxy. 

The surface brightness of an elliptical galaxy would make 
such a search difficult in the central regions. However, be- 
yond w 2 arcminutes, the surface brightness at V (in NGC 
4697, for example) is below the sky level (for HST) of 22-23 
Mag/arc -sec^ (Goudfrooij et al. 1984). In that limit, and us- 
ing the measured Chandra position of sources at such locations, 
such a wide (10 day) binary would be detectable at B or V in 
a 5 X 10"* s observation with the Advanced Camera for Surveys 
on HST. Of course, such optical counterparts are easy to find in 
the nearby bulge of M31, where V « 22. 

3. DISK INSTABILITY, TRANSIENT FRACTION, AND RECURRENCE 

TIMES 

If the disk is thermally stable and transfers matter to the NS at 
the supplied rate, there is nearly a one-to-one relation between 
Pj and the X-ray luminosity (see Figure 1) that would allow 
the orbital period distribution to be inferred from the luminos- 
ity distribution (Webbink et al. 1983). However, this argument 



fails because these wide binaries are almost always transient 
accretors. 

The usual condition used for the disk instability is that some 
portion of the disk is below the hydrogen ionization temper- 
ature (« 6500 K). Since the effective temperature of the disk 
decreases with increasing radius, it is usually adequate to check 
that the outer radius of the disk is below the critical temperature. 
In X-ray binaries, however, this temperature is fixed by irradia- 
tion from the central accreting source (van Paradijs 1996). The 
surface temperature is then estimated as (Dubus et al. 2001) 



CrjMic^ 
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where R is the disk radius, cr is the Stefan-Boltzmann con- 
stant, ?7 = 0.1 is the conversion efficiency for accretion onto 
the NS, and C measures the fraction of the X-ray luminosity 
that reaches the disk (and contains information on the irradi- 
ation geometry). We take C = 5.7 x lO""' (slightly larger than 
the preferred value, C = 5 x 10"^ of Dubus et al. 2001) so as to 
maximize the persistent fraction. The disk is assumed to extend 
out to 70% of the Roche lobe around the NS (King et al. 1997). 

Setting Tin- = 6500 K allows us to solve for the critical time, 
accretion rate, and orbital period at which the instability sets in 
as a function of the donor's mass. The critical mass transfer rate 
is then 
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which is used for the instability boundary shown in Figures 1 
and 2 as the dotted line. The critical orbital period is 



i'orb.cnt = 5.5 days Mj/" (2.3 -M2)- 
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and this is used for the instability boundary in Figure 2. The 
fraction of time spent transient is then found (numerically to 
better than a percent) as a function of the initial orbital period. 

Transient Fraction = -0.0500x^ -I- 0.2385x -I- 0.7655, (6) 

where x = logio (^//d)- For P,- in excess of a day, the binary 
is transient at least 75% of the time, and this does not include 
binaries containing black holes which are transient almost all 
the time (e.g. King et al. 1996). 

We estimate a maximum recurrence time for outbursts by 
assuming the outside-in model. We assume that the binary is 
quiescent until a critical surface density in the outer accretion 
disk (at radius R) has been reached (Dubus, Hameury,& Lasota 
2001), 
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where we assume that the viscosity parameter a = 0.1 and that 
there is no irradiation of the disk in quiescence. This expression 
gives a maximum recurrence time 
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that should be viewed with some caution given the uncertain- 
ties in the accretion disk physics. This recurrence time is 
100-10,000 years (depending on the initial orbital period; see 
Figure 2) and is consistent with current estimates from observa- 
tions of galactic ms radio pulsar/He white dwarf binaries (Ritter 
& King 2002). Our estimate is maximal because it is possible 
(and maybe even likely) that the LMXB may outburst before 
the outer disk reaches the critical mass, such as an outburst ini- 
tiated in the inner disk. 

4. COMPARISON WITH GALACTIC SOURCES 

One difficulty with comparing these evolutionary models 
with observations is the lack of long orbital period LMXBs in 
our galaxy. Many persistent sources have low enough orbital 
periods (^ 1-20 hours) that they are probably not accurately 
described by the evolutionary model used above. However, a 
lack of persistent sources at Po^b > 1 day is in agreement with 
our calculations and generally expected from previous work 
(for example King, Frank, Kolb & Ritter 1997). Figure 1 shows 
that such NS binaries spend most of their lifetime as transient 
rather than persistent sources (> 75% of the time). Though all 
these relations are consistent, it is still hard to claim that the 
instability criterion has been fully scrutinized by the galactic 
population. 

Indeed, the outburst durations of these systems might be so 
long that they appear to us as "persistent" sources (with ac- 
cretion rates most likely in excess of the mass transfer rate 
expected from binary evolution) until they suddenly disap- 
pear Indeed, one such source (and maybe others; Wijnands 
2002) has now been clearly identified: KS 1731-260 (Wij- 
nands et al. 2001; Rutledge et al. 2002). It was accreting 
at M « 3 X lO~^M0yr"' for at least a decade prior to its sud- 
den decline in luminosity by a factor of > 10^ (Rutledge et al. 
2002). Analysis of the immediate quiescent emission in terms 
of thermal emission from the NS surface yielded constraints 
on the recurrence time from 300-1,000 years (Wijnands et al. 
2001; Rutledge et al. 2002). Until this system's orbital period 
is found (current knowledge of the infrared counterpart points 
to a possibly evolved companion; Orosz, Bailyn & Whitman 
2001) we will not be able to say whether it is an LMXB of the 
type we are discussing. However, its existence does point to a 
likelihood of finding such similar systems in elliptical galaxies. 

One galactic binary that does have a long orbital period is 
GRO J1744-28 (11.8 days; Finger et al. 1996). This system 
is transient and was shown by Bildsten and Brown (1996) to 
have a companion with mass M2 « 0.334Mq and core mass 
Mc « O.22M0 (if the companion is a population I metallicity 
star). When this binary's evolutionary path is predicted and 



plotted, it is found to be transient (as it is) with a maximum 
outburst recurrence time of ^1,000 yrs (see dot-dashed line in 
Figure 2), consistent with a single outburst. 

5. CONCLUSIONS AND DISCUSSION 

We have shown that the expected mass transfer scenario for 
low mass X-ray binaries in the field of an elliptical galaxy 
(where the donors must be < Mq) that yield M « lO"^M0yr"' 
(as observed by Chandra) have thermally unstable disks. Such 
binaries are expected to be transient accretors with maximum 
recurrence times of 100-10,000 years for at least 75% of their 
lifetimes. When active, it appears likely that they would easily 
reach the NS's Eddington limit. 

Our inability to calculate the duration of the bright outburst- 
ing phase makes it difficult to predict how many such transients 
will be identified in repeated Chandra visits. If the duty cycle 
is bright for ten years out of a thousand (like implied for KS 
1731-260; Rutledge et al. 2002), and there are 100 such bina- 
ries, then a year later, roughly 10 would have faded into quies- 
cence to be replaced by new transients. This is not that far off 
from what was seen over a five month interval in Cen A (Kraft 
et al. 2001), where four sources (out of 120 exposed in both 
epochs) were either non-detectable (L^ < (1 -2) x 10^^ erg s"') 
or appeared bright (L^ > 4 x lO^^erg s"') in either epoch. There 
were a larger number (« 35) which showed significant variabil- 
ity (Kraft etal. 2001). 

Such comparisons are also complicated by the discovery that 
many of the X-ray sources in ellipticals are coincident with GCs 
(Sarazin et al. 2001; Angelini et al. 2001; White, Sarazin & 
Kulkarni 2002). Other possible accretion scenarios then include 
accretion from lower-mass main sequence stars that are formed 
in this dense stellar environment (see Rasio, Pfahl & Rappa- 
port 2000 for a recent discussion). These tight systems would 
most likely be persistent accretors, thus apparently reducing the 
fraction of transients amongst the detected sources. This might 
just be an impression however, as if the duty cycle of the field 
transients is really 1 %, then the underlying population is indeed 
larger than that in GCs. Namely, it is possible that the GC pop- 
ulation is only apparently a large fraction because LMXBs in 
GCs are persistently bright. The census of the much more nu- 
merous population of field transients can only be taken over a 
few thousand years. 

We thank G. Dubus, J. Irwin, H. Ritter, and C. Sarazin for 
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the Research Corporation. 
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Fig. 1 . — Mass transfer rate as a function of time for LMXBs with red giant donors. The solid lines are the evolutionary path of the binary given an initial orbital 
period when Roche lobe overfilling begins (evolution is from left to right). The dashed line is the boundary where the Hydrogen envelope of the donor has vanished, 
thus halting the evolution. The dotted line is the boundary at which the accretion disk instability begins. 




Fig. 2. — Estimated maximum recurrence time as a function of tlie orbital period. Tlie curves are tlie same as in Figure 1. The additional dot-dashed Une is 
predicted evolutionary path of GRO J1744-28. The filled circle denotes this binary's current orbital period of 1 1.8 days. 



